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Abstract 



The effects that the structure of a neutron star would have on the gravita- 
tional emission of a binary system are studied in a perturbative regime, and 
in the frequency domain. Assuming that a neutron star is perturbed by a 
point mass moving on a close, circular orbit, we solve the equations of stellar 
perturbations in general relativity to evaluate the energy lost by the system 
in gravitational waves. We compare the energy output obtained for different 
stellar models with that found by assuming that the perturbed object is a 
black hole with the same mass, and we discuss the role played by the excita- 
tion of the stellar modes. Our results indicate that the stellar structure begins 
to affect the emitted power when the orbital velocity is u ~ 0.2 c [vgw ~ 185 Hz 
for a binary system composed of two IAMq neutron stars). We show that the 
differences between different stellar models and a black hole are due mainly 
to the excitation of the quasinormal modes of the star. Finally, we discuss to 
what extent and up to which distance the perturbative approach can be used 
to describe the interaction of a star and a pointlike massive body. 

PACS numbers: 04.30.-w, 04.40.Dg 
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I. INTRODUCTION 



In a recent paper (to be referred to hereafter as Paper I) we have studied the grav- 
itational emission of a binary system by using the following perturbative approach: one of 
the two stars is assumed to be an extended body, whose equilibrium structure is described 
by an exact solution of the relativistic equations of hydrostatic equilibrium; the second star 
is a pointlike mass which induces a perturbation on the gravitational field and on the ther- 
modynamical structure of the extended companion. We modeled the extended star using a 
polytropic equation of state, we chose the parameters so that the radius and the mass were 
those of a plausible neutron star (NS), and we solved the equations of stellar perturbations 
to compute the power emitted in gravitational waves when the point mass moves on orbits 
of arbitrary eccentricity. 

In this paper we use the same perturbative approach to compute the power radiated 
in gravitational waves by different stellar models. The purpose of this investigation is to 
answer a number of questions that arise in studying the signals emitted by NS-NS binary 
systems during the last phases before coalescence. The first is about the role played by 
the internal structure of the star: how does the gravitational emission depend on the mass 
and radius of the star (and consequently on the equation of state (EOS) of dense matter)? 
In particular, what is the difference if the perturbed object is a black hole (BH)? The 
second is related to the possibihty of exciting the quasi normal modes of the star. Are 
the differences between stars and a black hole (if any) due exclusively to resonant modes 
or to other orbital effects? The third question concerns the domain of applicability of the 
perturbative approach: to what extent are the results of this approach applicable to describe 
the evolution of a true binary system, composed of two neutron stars of comparable mass? 
Are we still in a perturbative regime when the two neutron stars are, say, 3 or 4 stellar radii 
apart? 

The plan of the paper is the following. In Sec. II we shortly review the procedure 
we use to find the solution of the equations of stellar perturbations inside and outside the 
perturbed star, along the lines of Paper I. Since we are mainly interested in the last phases 
of the orbital evolution, we shall consider the orbit of the point mass as already circularized 
0. In Sec. Ill we discuss in detail the results of the integration of the perturbed equations 
for the different stellar models and for a black hole, both excited by the same process. In 
Sec. IV we examine the domain of applicability of the perturbative approach applied to 
binary systems, and in Sec. V we draw the conclusions. 

II. THE PERTURBED EQUATIONS 

In order to compute the radiation emitted by a a neutron star perturbed by a mas- 
sive point particle, we need to integrate the equations of stellar perturbations inside and 
outside the star, and this can be accomplished by using different gauges and formalisms. 
Inside the star we choose the Chandrasekhar- Ferrari gauge 0, which allows to decouple the 
equations for the gravitational perturbations from those describing the perturbations of the 
fluid. After expanding the perturbed metric tensor in tensorial spherical harmonics, choosing 
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ihis gauge the radial part of the polar metric perturbations is described by four functions, 
Nim{^^, r), Lim{uj, r), Vim{uj, r),Tim{uJ, r) , whereas only two functions, h'}^{uj, r), hj^{uj, r) , 
are needed to describe the axial part. As in Paper I, we integrate the equations of stellar 
perturbations in the frequency domain (cf. Paper I, Eqs. (2.2) and (2.4)) up to the surface 
of the star, where we construct the ZeriUi and the Regge- Wheeler functions, Zi^{u, r) and 
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(2.1) 



and their first derivatives with respect to r. The integration of the equations of stellar 
perturbations in the interior (together with the TOV equation) is done using an adaptive 
step Runge-Kutta method. The relative tolerance of the integrator is set to 10~®. We use 
as integration variable the logarithm of the pressure, which allows us to approach the stellar 
surface gradually, and to determine the radius with the required accuracy. 

The two functions in (|2.1|) are needed to compute the radial part, '^im{uj,r), of the 
perturbation of the Weyl scalar which is defined as 
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Outside the star we integrate the inhomogeneous Bardeen-Press-Teukolsky BPT equation 
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where A = — 2Mr, and the source term Ti^iuJ, r) is that appropriate to describe the 
point mass mg moving on a given orbit around the star. In Paper I we discussed how to 



construct the solution of eq. (|2.4| ) in the general case of elliptic orbits. In this paper, since 
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we are interested mainly in the last phases of the evolution of binary systems, when the 
orbit has already been circularized, we focus on circular orbits and give the explicit solution 
of the equations in that case. If the mass tuq moves on an orbit of radius Rq the geodesic 
equations give 
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where the dot indicates differentiation with respect to proper time. E is the energy of the 
particle per unit mass, and uo^ is the keplerian orbital frequency 
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The source term can be written as 
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The functions sUim are given by 

-2'KJnin + l) ^ 



0'S'itn( — , 0) oUirn +-1 'S'f„( — , 0) _lf//m + -2'S'/*f„( — , 0) ^2^^/ 



(2.6) 



(2.7) 



oUlm 
-lUlm 



5(r - i?o) 



6{r - Ro) (27T V2^ m 7 Aq i?o) - S'{r - Rq)— [2m V2^ 7 ujk RI 



-2Uim = S{r - Ro) 
+S'ir - Ro 



IS? 
.3 - ^0 



ivr m uoj^ 7 



Rl 



r 



TT c<j^ 7 Ro 



+ 



Ar^ (27ii rauo]^ 7 Rq) + A^ [Ah uj\ 7 Ro 



+ 



+5'\r - Ro) [aV (vr col 7 R, 



where the prime indicates differentiation with respect to r, and a subscript means eval- 
uation at r = Rq. As in Paper I, the solution of Eq. ( |2.4| ) which satisfies the boundary 
conditions of pure outgoing radiation at radial infinity, and matches continuously with the 
interior solution can be found by the Green's function technique, and the amplitude of the 
wave emerging at radial infinity can be shown to be (cf. Paper I, Eqs. (4.3)-(4.7)) 
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where Wim{^^) is the Wronskian of the two independent solutions of the homogeneous BPT 
equation 
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and "^i^ and \E';^ satisfy the following equations 



]LBPT^,°(a;,r) = 



'lLBPT^,^(^,r) = 
^iliuJ.R) = ^imiuj,R) 



1 ', 

Im 



(2.10) 



00, R) 



where ^bpt is the differential operator on the left hand side of the BPT equation, and 
^im{uj,R) is the function constructed according to Eq. (^73|) , in terms of the perturbed 
functions computed in the interior, and evaluated at the surface of the star. The integral in 
Eq. (|2.8|) can explicitly be evaluated by making a repeated use of the following property of 
the 5-function 
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If we define 

the time-averaged energy-flux 
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can be written in terms of 
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In order to evaluate and the BPT equation is integrated with the same adap- 
tive Runge-Kutta used for the equations in the interior. We must remark that close to a 
resonance the solutions need to be computed very accurately, since the Wronskian is the 
difference between two terms that almost cancel each other. When required, the tolerance 
in the integration is decreased until convergence is reached. Typically, computing the power 
emitted at a given orbit takes about 10 s on a PC with a Pentium IV IGHz processor. 

Since the orbital frequency is related to the orbital velocity v and to the semilatus rectum 
p = Ro/M hj the following relations 

V = {MuKf'^ = ^, (2.15) 

the energy flux E^ can also be considered as a function of v or p. 

In the following, we shall normalize E^ to the Newtonian quadrupole energy flux 
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III. COMPARING THE GRAVITATIONAL FLUX OF STARS AND BLACK 

HOLES 



In our analysis, we consider five models of stars with a polytropic EOS, using two values 
of the polytropic index, n = 1 and n = 1.5. The values of the central density, the 
ratio ao between the energy density and the pressure at the centre, the mass, the radius 
of the star, and the compactness M/ R (in geometrical units) are given in Table |. We have 
chosen polytropic models for simplicity, but without loss of generality, because the differences 
in the gravitational flux are expected to depend more on global properties such as mass, 
radius, average density, or compactness, rather than on the specific matter distribution. The 
parameters we choose encompass a reasonable range of stellar models (radius ranging from 
9 to 15 km), and the polytropic exponents, r = l + l/n = 5/3 and T = 2, cover most 
of the range of structural properties obtained with realistic EOS's. 

A. Power radiated and resonant modes 

For each model, we integrate the equations of stellar perturbations as described in Sec. 
II, assuming that a point like mass, mo, is moving on a circular orbit of radius Rq with orbital 
velocity v, and we compute the energy flux emitted in gravitational waves normalized to 
the Newtonian quadrupole energy flux, P = E^/E^, where E^ and E'^ are given 
in Eqs. ( |2.14| ) and ( |2.16| ). In Ref. P^TU] the gravitational emission of a Schwarzschild black 
hole perturbed by a massive particle in circular orbit was studied numerically in great detail. 
Our results for the black hole agree with them with an accuracy of (at least) one part in 
10®. In order to allow an easy comparison, and for future reference, we show in Table y the 
contribution of the 2 < / < 5 multipoles to the total gravitational power. The results for 
the black hole can directly be compared to Table II in Ref. p!0| . 

In Fig. 1^ we plot the normalized energy flux, P{v), as a function of the orbital velocity, 
for the models of star we have considered, and for the black hole. P{v) has been obtained 
by adding the contributions of different Us and m's, with 2 < / < 7. As discussed in paper 
I, a multipole of order / contributes to the total power as a correction of order p'^~K By 
truncating our multipole expansion at / = 7 we are incurring in a relative error of order 
= v^"^. It should be mentioned that P{v) is independent of the mass mo. For the black 
hole and for the stellar model D, the plots extend up to the velocity v = 0.408, which 
corresponds to the Innermost Stable Circular Orbit (ISCO), -Ro = 6M; for model A, 
B and C the ISCO would be in the interior of the star, therefore the plots are truncated 
at a smaller velocity, which corresponds to mg reaching the stellar surface. Sharp peaks 
appear if the central object is a star: they correspond to the excitation of the fundamental 
quasi-normal modes of the star for different values of the harmonic index /. In the case of 
model B the first p-mode for Z = 2 is also excited. In Table |ITT| we show the values of the 
radius Rq, of the dimensionless orbital velocity v, and of the keplerian frequency of the 
orbit that corresponds to the excitation of the fundamental modes of the star for different Us 
for the considered stellar models. The corresponding frequencies of the /-mode are given in 
the last column. From the analytical form of the stress-energy tensor ( p.7|) , it is easy to see 
that, for each assigned /, a mode of the star is excited when the orbital frequency satisfies 
the resonant condition 
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(3.1) 



where i>i is the mode frequency. Table III shows that the frequency of the f-mode increases 
with /; however, eq. ( |3.1| ) shows that, for instance, the orbital velocity that corresponds 
to the excitation of the /-mode for / = 3 is lower than that needed to excite the /-mode 
for / = 2; this means that in the process of coalescence of the "binary system" formed by 
the star and the point mass mo, the /-mode for / = 3 is excited before the / = 2 /-mode, 
and similarly the / = 4 f-mode is excited before that for I = 3, while the gravitational wave 
frequency is higher. The peaks corresponding to higher / are narrower, so that peaks for 
/ > 4 are difficult to locate, even with the use of an extremely refined grid. From newtonian 
theory we know that the /-mode frequency scales with the square root of the average density 



of the star. This is true also in general relativity and the dependency on J M/E? is still 



linear [|T^]. This explains why for a chosen value of the polytropic index n the peaks for 
more compact stars occur at higher frequency, i.e. at higher f , (for instance, compare in 
Fig. m the two curves for model B and D, for which n = 1 and i? = 15 km and R = 9.8 
km, respectively). Since the peaks corresponding to the mode excitation are very high, the 
scale chosen on the vertical axis of Fig. |l| makes the response of the black hole to appear 
as a fiat line. The reason is that, since the frequency of the lowest quasi-normal modes of a 
black hole are higher than those of a star with the same mass, the circular orbit that would 
excite them would have a radius smaller than QMbh- Thus, in the range of v considered in 
Fig. I the energy fiux emitted by the black hole is due essentially to the orbital motion. In 
Fig. 0, we show a zoom of Fig. |I| restricted to the region v < 0.28, which is far enough from 
the resonant orbits (except that for model A). In this case we can appreciate the differences 
between the emission of different stellar models and that of a black hole. If the orbital 
velocity is smaller than 0.16 all curves are practically indistinguishable. 

Fig. ^ shows that the normalized energy fluxes emitted by different stellar models have 
a different slope, and are always larger than the flux emitted by the black hole. The curve 
for model E (n = 1.5, -Rq = 9 km) is practically indistinguishable from the black hole curve; 
that for model D (n = 1, i?o = 9.8 km) is also very close to the black hole result. At flrst 
sight we may relate this behavior to the fact that these stars are more compact than models 
A,B,C; however, the steepest raise of the curves of models A,B,C could also be a marginal 
effect of the resonances, or may be due to a different coupling between the orbital motion 
and the stellar structure. 



B. A harmonic oscillator model 

In order to better understand the underlying physical picture, we shall use a toy model 
which has been employed in many contexts to study the effects of stellar resonances (see for 
instance Ref. ||12[)- Since a star oscillating in a quasi-normal mode of complex frequency 
Uo — iuJi emits a gravitational wave of amplitude ^ Ce~^(^o - i^^i)t ^ jt 

can be modeled as 

a harmonic oscillator which satisfles the differential equation X + 2uiX + {uJq + ujf)X = . 
Here the amplitude X is assumed to be normalized to some reference amplitude, for instance 
to the amplitude of the newtonian quadrupole, ^4^, which we deflne in terms of the newtonian 
energy flux of Eq. ( p.l6|) as follows 
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E- = ^^\A-[u)\\ (3.2) 

The particle orbiting the star acts as a driving force proportional to cj^e"*'^*, and the ampli- 
tude of the forced oscillation is found by solving the inhomogeneous equation 

X + 2u:iX + [ul + uDX = buj^e-'''' ; (3.3) 

the solution can be written in the form X{t) = X{uj)e^^'^^ , where 

= 2~K^o- • (3.4) 

We shall now assume that, near a resonance, the total amplitude of the wave emitted by the 
perturbed star is the sum of two contributions, one due to the orbital motion, the other to 
the star pulsating in its quasi-normal mode, i.e. 

AU^) = A''{uj)[l + X{u)]. (3.5) 

According to this model, the normalized energy flux P^^s can be written as 

Presij^) has a maximum in ujq and a minimum in ° ; the height of the maximum is 

Vl — 6 

- — 5^ and that of the minimum is — \ „ (up to corrections of order OiiuJi/ ujqY)). 

For uj « uq, Presi^) tcuds to 1, while after the resonance it tends to (1 — 6). In order to 
check how good is the harmonic oscillator model to describe the behavior of the star near a 
resonance, the quantity Pres(i^) has to be compared with the normalized energy flux, P{uj), 
computed by integrating the equations of stellar perturbations. In Fig. ^ we plot Pres(t^) 
(continuous line) and the numerically evaluated Pi=2,m=2{'-^) (diamonds) as a function of 
ujM, in a frequency region near the fundamental mode for / = 2, m = 2, for the stellar 
model C. The parameter b is found from the locations of the maximum and the minimum, 
i.e. 6=1 — {wmax/^miny, and Ui is chosen in such a way that the amplitude of the maxima 
of the two functions coincide. The results are in very good agreement (relative differences 
are of order 10~^) except very close to the minimum, where the function Pi=2,m=2{^) goes 
to zero. Incidentally, it should be mentioned that the evaluation of the exact value of the 
minimum for each / is not so important, as far as the total energy flux is concerned; indeed, 
the total flux P{uj), evaluated at the frequency which corresponds to a minimum for a given 
/, has a finite value, because it is the sum of contributions from different Vs and m's, and it 
is dominated by the other multipoles. 

The simple analytical model works surprisingly well in the whole range of each resonance, 
capturing the parabolic behavior in the region |a; — cuop <C 1, used by several authors to 
calculate damping times, and also describing the global behavior of the normalized flux with 
high accuracy. We found that the same procedure can satisfactorily be applied to the other 
stellar models, and to resonances corresponding to higher values of I. The values of the 
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frequencies where the maximum and the minimum occur, the forcing amphtude b and the 
imaginary part of the frequency of the fundamental mode, Ui, are summarized in Table 
for the five stellar models. Note that the cUjM scales approximately as 10~^', and b as 
10-'. 

The agreement between the numerical result and the toy model suggests the following 
question: can the effects of a resonance be subtracted from the global signal ? The answer 
is shown in Fig. ^ where we plot the (/ = 2, m = 2) contribution to the flux, P(f )i=2,m=2, 
with a solid line, and the result of subtracting the effect of the resonance with dashed 
lines, versus the orbital velocity. The resonance-free curve has been obtained dividing the 
normalized energy flux emitted by the perturbed star by the model given in (|3.6| ), i.e. 
P{v)i=2,m=2l Presiv). Now the questiou of whether or not the different slope of the curves 
plotted in Fig. |^ can be attributed to the effect of resonant excitation of a stellar quasi- 
normal mode can be answered positively. The increase in the energy output at orbital 
velocities of about v = 0.2 is just an effect of the resonance tail. However, it is not the 
only effect. Indeed, since we are normalizing the emitted flux to the newtonian quadrupole 
flux, if the signal emitted by the perturbed star would be, as we assume in eq. (p.5[), the 
sum of a term due to the orbital motion (the quadrupole) and a term due to the resonance, 
the dashed line in Fig. ^ should be a horizontal straight line. Conversely, it is a slightly 
decreasing function of v. In order to see if this is a general feature, we have also calculated 
P{v)i=2,m=2/ Pres{v) foY the stcllar models A, B and D. The results are plotted in Fig. |^. 
For comparison, we also plot -P(f )z=2,m=2 for a black hole and for the stellar model E, 
computed without subtracting the contribution of the resonances. The reason is that we 
cannot subtract this contribution because, as mentioned before, black hole quasi-normal 
modes cannot be excited before the ISCO and similarly, to excite the f-mode of model E 
the point mass should move on an orbit with radius smaller than the stellar radius. We see 
that after the subtraction there is still a difference between different EOS's; in particular, 
if we classify the models A,B,C,D according to their stiffness (where a stiffness indicator is, 
for example, the speed of sound at a given density), we find that the stiffer the EOS is, the 
smaller is the orbital contribution to the total emission. We must remark, however, that 
when the effect of the resonances is included, this trend is inverted; stiffer models emit more 
energy, because their resonant frequencies are smaller and their effect is more pronounced. 

The influence that the structural effects discussed above may have on the detection of 
these signals by terrestrial detectors, should be estimated by computing the phase evolution 
of the signal; this can be done by taking into account the evolution of the orbit due to 
radiation reaction effects. Although this problem will be specifically addressed in a separate 
paper, we can anticipate that, depending on the stellar model, most of the phase shift is 
accumulated during the last 50 cycles before merging; this effect is not likely to be detectable 
by the first generation of LIGO and VIRGO detectors, unless we are so lucky that a NS- 
NS coalescence occurs closer than 10 Mpc. However, the situation would be different if 
a detector like EURO |T^, recently proposed, would be constructed, because it would be 
very sensitive in the frequency range (600-1200) Hz where the structural effects become 
significant. Quantifying precisely the differences between different signals and templates is 
out of the scope of this paper, and requires a careful calculation of the overlaps of the "true 
signal" with those that are currently used as templates. This work is in progress and will 
be reported in a forthcoming paper. 
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IV. ON THE VALIDITY OF THE PERTURBATIVE APPROACH 



The theory of perturbations of black holes excited by a point particle developed since the 
early seventies was based on the assumption that the mass of the particle is much smaller 
than the black hole mass; under this condition, the effect of the particle on the black hole 
is that of inducing a small perturbation on the equilibrium configuration, and the stress 
energy tensor of the mass tjiq can be considered as a source for the perturbed equations. In 
addition, since mg does not affect the background geometry, it will move on a geodesic of 
the unperturbed spacetime. The same assumption can be done to study the perturbations 
of a star, as we do in this paper. However, since one of the purposes of our study is to have 
an insight on the last phases of the coalescence of a true binary system composed of stars of 
comparable mass, we may ask the following question: how big can the mass mo be, in order 
the deviation it induces on the gravitational field and on the thermodynamical structure 
of the companion star to be considered as a perturbation? Or, we can formulate the same 
question in a different way: given a mass mg, not necessarily small, moving on a circular 
orbit around a star, up to which distance its effect on the companion can be considered as 
a perturbation? Q To answer this question in Fig. ^ we plot the fluid and the gravitational 
perturbed functions inside the star, for / = 2,m = 2 (the most significant contribution) and 
for the stellar model D, assuming that the point mass is moving on a circular orbit at a 
distance Rq = 3R from the star. In the upper panel we plot the radial component of the 
Lagrangian displacement normalized to the stellar radius and the Lagrangian perturbations 
of the density and of the pressure normalized to their equilibrium values. In the lower panel 
we plot the perturbations of the metric 2N[r) and 2L(r), that have to be compared to unity 
(cf. Paper I, Eq. (2.1)). All quantities are normalized to the ratio rrio/M. Similar results 
are obtained for the other stellar models. Fig. |] shows that, even if we assume that the 
two masses are equal, all perturbations belonging to the fluid or to the gravitational field 
are small compared to the corresponding unperturbed functions. Thus, the perturbative 
approach holds even when the two bodies have comparable mass and are as close as three 
stellar radii. Corrections of the same order of magnitude have been found in full nonlinear 
stationary solutions of binary neutron stars M 



V. CONCLUDING REMARKS 

In this paper we have studied in a full general relativistic approach how the internal 
structure of a neutron star affects its gravitational emission when the star is perturbed by a 
close, orbiting companion. Our study has been done by using a perturbative approach in the 
frequency domain, and the massive object which perturbs the star has been assumed to be 
a point particle of mass mo. We have considered five models of star with global properties 
- like compactness, average density, stiffness of the EOS - that encompass a wide range of 



We neglect the fact that if rriQ is large it will not move on a geodesic of the unperturbed 
spacetime. This point will be discussed in the concluding remarks. 
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stellar properties, and we have compared their behavior with that of a black hole excited by 
the same source. 

The results we find can be summarized as follows. The internal structure of the star 
affects the emitted fiux of gravitational waves only when the orbital velocity of rriQ is (ap- 
proximately) higher than one fifth of the speed of light, i.e. when mo is very close to the 
central star. If the mass of the star and mo were comparable and about 1.4 Mq, this velocity 
would correspond to an emission frequency ~ 185 Hz. For lower values of v (larger orbital 
separation) neutron stars and black holes behave in the same manner and the gravitational 
fiuxes they emit are practically indistinguishable. 

The difference in the energy flux emitted by a star and a black hole arises mainly because 
of the excitation of the fundamental mode of the star for different /'s, which can be modeled 
extremely well in terms of a suitably deflned harmonic oscillator. The results for model A 
suggest that the emission properties of low mass neutron star binary systems would deviate 
signiflcantly from those of black hole-black hole binary systems. 

Once the effect of the resonant excitation of the quasi-normal modes of the star are 
subtracted, we flnd that there is still a residual difference between the emission of neutron 
stars and black holes, and that stars with stiffer EOS emit less energy than a black hole. 
However, this trend is inverted when the effect of resonances is included. 

It should be stressed that the region where corrections due to the EOS begin to be 
signiflcant is the same region where high order Post Newtonian (PN) corrections play a 
signiflcant role. Thus these effects should be included in constructing templates to be used 
in the data analysis of gravitational detectors. In addition, the corrections introduced by 
the tail of the resonant f-mode is comparable, if not larger, to that introduced by the high 



order PN corrections [15 



The problem of mode excitation during the latest phases of the evolution of a binary 
system, and the consequences it has on the emitted gravitational signal, have previously been 
studied by Ho & Lai []TB[ by a newtonian approach. In particular, as far as the excitation of 
the f-mode is concerned, they consider the case of a rotating, incompressible model of star, 
which is very close to our model E, whereas we use a relativistic perturbative approach and 
focus on the differences which may arise because of different EOS's. Apart from this, we 
basically reach similar conclusions about the order of magnitude of the resonant effects and 
on the regime where they may play a role. 

It is interesting to discuss to what extent can we extrapolate our results to the case 
when the mass mg is comparable to that of the central star, i.e, to simulate a true coalescing 
binary system. Indeed, the results of Sec. IV indicate that the perturbative approach holds 
also when two bodies of comparable mass are very close, even closer than 3 stellar radii. A 
naive generalization of the energy fluxes we obtain could be that of rescaling the amplitude 
of the energy flux by the appropriate value of mo, and to change the orbital frequency 



uk = yM/Rl (and consequently the emission frequency u>gw = itlujk) by replacing the 
mass of the star M with the total mass of the system Mf = M + mo; it is worth mentioning 
that, in this case, given a certain orbital radius Rq, the orbital frequency corresponding to 
that radius would be higher than that of the point particle on the same orbit. This can easily 
be done, but the results should be considered only as an indication of what may happen in 
reality. Indeed, in order to correctly generalize the results, the perturbative approach should 
be improved in many ways: the geodesic equation which are used to describe the motion of 
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mo should be replaced by the correct equations of motion of a two-body system, as seen in 
the coordinate system centered on one star; these equations are now known at the 3.5 PN 
level |]T^, and they would change the stress energy tensor we put on the right-hand side of 



our equations; the problem of chosing the reference frame in which to compute the radiated 
power as seen by a distant observer is a very delicate one, and in order to compare the results 
of the perturbative approach for equal masses with those of the post-newtonian approach one 
should change to a frame centered in the center of mass of the binary system; in computing 
the orbital evolution due to radiation reaction effects it is customary to use the adiabatic 
approximation which assumes that the timescale of the orbital evolution is larger than the 
orbital period; we think that the role of this assumption has to be investigated very carefully 
during the last few cycles before coalescence. All of these problems are important: we are 
looking for very small effects which arise just before coalescence, and they may introduce 
corrections of a few percents in the emitted power, which are of the same order as the effects 
of stellar structure we discuss in this paper. Finally, one should certainly consider rotating 
stars. Rotation may change the situation in a positive direction, because we know that it 
has the effect of lowering some of the mode frequencies, and this may enhance the marginal 
effect of the mode excitation which we see in the tail of the signals emitted by compact stars. 

In future papers we plan to implement all these effects in our perturbative approach, 
with the final goal of providing accurate templates which may be used in the data analysis 
of future high sensitive detectors operating in the kHz region. 
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TABLES 



TABLE L Parameters of the polytropic stars we consider in our analysis: the polytropic index 

n, the central density, the ratio ao = cq/pq of central energy density to central pressure, the mass, 
the radius and the ratio M/R (ag and M/R are in geometric units). The central energy density 
is chosen in such a way that the stellar mass is equal to IAMq, except for model A, the mass of 
which is about one solar mass. 



Model number 


n 


Pc (g/cm^) 


"0 


M (Mo) 


R (km) 


M/R 


A 


1.5 


1.00 X 10^^ 


13.552 


0.945 


14.07 


0.099 


B 


1 


6.584 X 10^^ 


9.669 


1.4 


15.00 


0.138 


C 


1.5 


1.260 X 10^5 


8.205 


1.4 


15.00 


0.138 


D 


1 


2.455 X 10^5 


4.490 


1.4 


9.80 


0.211 


E 


1.5 


8.156 X 10^5 


2.146 


1.4 


9.00 


0.230 



TABLE XL Decomposition of the gravitational luminosity {M/mo)'^E^ into multipole con- 
tributions for p = 10. The BH results can be compared directly to those of Table II in Ref. 
[10]. 





m 


BH 


A 




B 




C 




D 




E 




2 


2 


0.536879(-4) 


0.339696( 


-04) 


0.220321(- 


06) 


0.667680(- 


04) 


0.542028(- 


04) 


0.537526(- 


-04) 


2 


1 


0.193161(-6) 


0.194489( 


-06) 


0.193693(- 


06) 


0.193491(- 


06) 


0.193244(- 


06) 


0.193189(- 


-06) 


3 


3 


0.642607(-5) 


0.588917( 


-05) 


0.599186(- 


05) 


0.564656(- 


05) 


0.643160(- 


05) 


0.642658(- 


-05) 


3 


2 


0.479591 (-7) 


0.481018( 


-07) 


0.479889(- 


07) 


0.479742(- 


07) 


0.479605(- 


07) 


0.479594(- 


■07) 


3 


1 


0.571489(-9) 


0.712057( 


-09) 


0.587559(- 


09) 


0.577923(- 


09) 


0.571763(- 


09) 


0.571520(- 


■09) 


4 


4 


0.953958(-6) 


0.918437( 


-06) 


0.944728(- 


06) 


0.948295(- 


06) 


0.954066(- 


06) 


0.953966(- 


-06) 


4 


3 


0.877874(-8) 


0.879191( 


-08) 


0.878020(- 


08) 


0.877939(- 


08) 


0.877878(- 


08) 


0.877876(- 


-08) 


4 


2 


0.526223(-9) 


0.116026( 


-08) 


0.532315(- 


09) 


0.528287(- 


09) 


0.526252(- 


09) 


0.526227(- 


-09) 


4 


1 


0.145758(-12) 


0.145966( 


-12) 


0. 145781 (- 


12) 


0.145768(- 


12) 


0.145759(- 


12) 


0.145759(- 


■12) 


5 


5 


0.152415(-6) 


0.149988( 


-06) 


0.152100(- 


06) 


0.152269(- 


06) 


0.152418(- 


06) 


0.152415(- 


■06) 


5 


4 


0.149211(-8) 


0.149338( 


-08) 


0.149219(- 


08) 


0.149214(- 


08) 


0.149212(- 


08) 


0.149212(- 


-08) 


5 


3 


0.182910(-9) 


0.168856( 


-09) 


0.184547(- 


09) 


0.183313(- 


09) 


0.182911(- 


09) 


0.182910(- 


-09) 


5 


2 


0.381934(-12) 


0.382239( 


-12) 


0.381952(- 


12) 


0.381942(- 


12) 


0.381935(- 


12) 


0.381935(- 


-12) 


5 


1 


0.236763(-15) 


0.252126( 


-15) 


0.237281 (- 


15) 


0.236943(- 


15) 


0.236765(- 


15) 


0.236764(- 


■15) 
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TABLE III. In this table we give the values of the radius, orbital velocity and keplerian 
frequency {i'k) of the circular orbits which correspond to the excitation of the fundamental mode 
of the considered stars for the first relevant multipoles, whose frequency is given in the last column. 
For model E we do not give these data for I = 2 because in order to excite the corresponding mode 
Rq should be smaller than the ISCO (6M). 



Model 


/ 


Ro (km) 


V 


UK (Hz) 


Uf (Hz) 


A 


4 


31.8 


0.255 


567 


2260 




3 


29.0 


0.267 


651 


1953 




2 


26.0 


0.282 


767 


1534 


B 


3 


22.9 


0.300 


626 


1879 




2 


21.0 


0.313 


711 


1422 


C 


3 


21.2 


0.312 


702 


2105 




2 


18.9 


0.331 


835 


1671 


D 


3 


15.5 


0.365 


1119 


3358 




2 


14.1 


0.383 


1296 


2593 



E 3 13.5 0.391 1379 4138 



TABLE IV. 


Parameters of the resonances 


(p^) corresponding 


5 to the fundamental 


mode of the 


most relevant /-multipoles. 


All frequencies are 


normalized to the 


mass of the corresponding star. 


Model 


I 


{ujM)max 


{ujM)min 


b 


LOiM 


A 


4 


0.066218 


0.066239 


0.000645 


9.67x10-1° 




3 


0.057314 


0.057495 


0.006272 


1.32x10"'^ 




2 


0.044823 


0.046273 


0.061682 


1.81x10"^ 


B 


3 


0.081420 


0.081645 


0.005618 


1.93x10-"^ 




2 


0.061581 


0.063363 


0.055367 


1.90x10-5 


C 


3 


0.091210 


0.091402 


0.004208 


2.09x10"'^ 




2 


0.072299 


0.074036 


0.046373 


2.47x10-5 


D 


3 


0.145709 


0.145895 


0.002555 


6.17x10"'^ 




2 


0.112454 


0.114287 


0.031836 


5.77x10-5 


E 


3 


0.179606 


0.1796935 


0.000971 


4.91x10""^ 
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FIG. 1. The normalized energy flux, P{v), is plotted as a function of the orbital velocity for 
the stellar models given in Table | and for a black hole. For model D and for the black hole the 
curves extend up to the velocity which correspond to the ISCO, whereas for the other models they 
stop when the mass tuq reaches the surface of the star. The sharp peaks indicate that, for different 
values of the harmonic index /, the fundamental quasi- normal modes of the star are excited if the 
orbital frequency satisfies the resonant condition the curve of the stellar model B has a peak 

at high V which correspond to the excitation of the first p-mode for / = 2. The most compact model 
E is not shown in the figure because at this scale it is indistinguishable from the black hole. 
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0.16 0.18 0.20 0.22 0.24 0.26 0.28 

V 

FIG. 2. The normalized energy flux, P{v), is plotted as in Fig. |l], but for a smaller orbital 
velocity range, such that the peaks due to the excitation of the stellar modes are excluded. 
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0.070 0.072 0.074 0.076 

(jjM 

FIG. 3. The normalized energy flux emitted by the stellar model C, is computed by using 
the toy model (solid line) and by integrating the equations of stellar perturbations for / = 2 and 
m = 2 (diamonds), and plotted as a function of the dimensionless orbital frequency uM near the 
resonance of the /—mode. The agreement is excellent. 
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FIG. 4. The I = m = 2 contribution to the normahzed energy flux, Pi=2,m=2{v), is plotted 
as a function of the orbital velocity (solid line), and compared to the result of subtracting the 
contribution of the resonance (dashes) using the analytical model, as described in the text. The 
data refer to the stellar model C. 
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FIG. 5. In this figure we plot the normalized energy flux emitted by the four models of star for 
I = 2 and m = 2 versus the orbital velocity; the contribution of the resonances has been removed 
by dividing Pi=2,m=2{v) by Pres, as described in Sec. III. For comparison, we plot the normalized 
energy flux emitted by the black hole (BH) and by the stellar model E {E{res)); in this case we 
do not subtract the contribution of the resonant excitation of the quasi-normal modes because in 
order to excite these modes the point particle should move on an orbit of radius smaller than 6M 
(for the black hole), or smaller than the stellar radius (for model E). 
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FIG. 6. This figure refers to the stellar model D perturbed by a point mass moving on a 
circular orbit of radius Rq = 3R. All quantities are normalized to the ratio mo/Ad, and only the 
I = 2,m = 2 contribution is shown, which is the most significant. In the upper panel, the radial 
component of the lagrangian displacement normalized to the radius of the star, ^r{f)/R, and the 
lagrangian perturbations of the density and of the pressure normalized to their equilibrium values, 
Ap(r)/p{r) and Ae(r)/e(r), are plotted as a function of the radial distance from the center of the 
star. The metric perturbations 2N and 2L plotted in the lower panel, have to be compared to 
unity, since they are a measure of the deviations of the perturbed metric functions z^(r) and //2(?") 
with respect to their equilibrium values (cf. Paper I, Eq. (2.1)). 
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